Effective interactions between star polymers 
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We study numerically the efltective pair potential between two star polymers with equal arm 
lengths and equal number / of arms. The simulations were done for the soft core Domb- Joyce 
model on the simple cubic lattice, to minimize corrections to scaling and to allow for an unlimited 
number of arms. For the sampling, we used the pruned-enriched Rosenbluth method (PERM). We 
find that the potential is much less soft than claimed in previous papers, in particular for / ^ 1. 
While we verify the logarithmic divergence of V{r), with r being the distance between the two cores, 
predicted by Witten and Pincus, we find that the Mayer function for / > 20 is hardly distinguishable 
from that for a Gaussian potential. 
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Interactions between polymers in diluted solutions are 
of interest for several reasons, not the least because they 
influence both the equilibrium and the rheological prop- 
erties of complex fluids. In early work by Flory et al. jlj it 
was suggested that polymer coils can be approximated by 
hard spheres, but this was shown to be wrong in . Since 
then it is well understood that both linear and branched 
polymers are soft in the sense that they can penetrate 
each other, and that the effective potential is a rather 
smooth function of their distance. As shown in 0|, 
this can have dramatic effects on the phase diagram for 
semi-dilute solutions of star polymers, and can - with the 
effective potentials assumed by these authors - lead to a 
multitude of novel phases. 

When discussing effective potentials between polymers 
- be they linear or star-shaped - one has to distinguish 
between U (r) where r is the distance between the two 
centers of mass, and V{r) where r is the distance between 
the two central monomers. In both cases, the potential 
is defined by 



exp(-/3;7(r)), exp(-/3y(r)) = Z^'^\r) / [Z^^^ 



(1) 



where Z*-^) is the partition function of a single polymer, 
while Z*-^) (r) is the partition function of two polymers 
with fixed distance r. Finally, (3 = X/ksT is used to give 
V{r) the usual dimension of a potential, although any 
temperature is of course dummy for an a-thermal system 
as in the present case. For ease of writing, we shall set 
/3 = 1 in the following. Finally, all data shown in the 
following refer to lattice simulations with r = (r, 0,0), 
but we checked in a few cases that distances not parallel 
to one of the coordinate axes gave basically the same 
results. 

In the following we shall only discuss the case where 
the number / of arms is the same for both polymers (and 
might include the case / = 2 describing linear chains), 
and each arm has the same length N . Even if not spelled 
out explicitly, the main point of H HQ is that, for large 
/, the potential V{r) is more relevant than U{r) for equa- 
tions of state of semi-dilute or dense solutions, and that 
V{r) is very different from U{r): While the latter is essen- 
tially Gaussian, the former has a more complex structure 
with a Yukawa tail at large r. We will show in the fol- 
lowing that at least the second claim is not correct, and 
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FIG. 1: Logarithms of the effective potential U{r) with r 
being the center-of-mass distance, for star polymers with 
/ — 2,3,4,5,6,8, and 10 arms (bottom to top), plotted 
against (r/Rg)^ where Rg is the gyration radius of a single 
star. Arm lengths are N = 400 monomers. Linear curves 
would correspond to Gaussian U{r). 



that V{r) can also be approximated by a Gaussian for 
most practical purposes. 

The center-of-mass potential U (rHs well known to be 
approximately Gaussian for linear 0, 0, Q, polymers. 
For star polymers there are much fewer computations 
0, so we present in Fig. 1 our own results which clearly 
indicate that U{r) is roughly Gaussian, too. Notice that 
the deviation from a Gaussian at small r (i.e. the upward 
bending in Fig. 1) is practically irrelevant for > 8 arms 
per star, since it occurs only when e~^^'^'> < 10"^. 

The data in Fig. 1 as all data in this paper were ob- 
tained for the soft repulsion Domb-Joyce model ^3 
the 'magic' value v* = 0.6 of the repulsion parameter, on 
the simple cubic lattice. This model was chosen because 
it leads to minimal corrections to scaling and it allows 
an arbitrary number / of arms to be attached to a single 
central site 0, The simulations were made with the 
PERM algorithm 'IS'l, adapted for star polymers as de- 
scribed in (1L] . The partition sum Z'^^^(r) was estimated 
as usual (e.g. 0) by simulating two independent stars si- 
multaneously, and computing their overlaps at different 
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TABLE I: Main results. The numbers in brackets are single 
standard deviations in the last digit. ^/ is defined by » 
AfN^'^, bf is obtained through Eq. ^ from the data of [ill . 
?7(0) is the effective potential when the two centers of mass 
coincide, and a/ is defined in Eq. (HJ, c/ and df are defined 
in Eq. and Tf is defined in Eq. IIOII . We do not quote 
errors for the latter four, since they are strongly correlated 
and individual error estimates would not make sense. 
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(7(0) 
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c/ 


df 
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0.2902(2) 


0.815(2) 


1.791(2) 


1.869 


0.372 


0.405 


4.5 
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0.3587(2) 


1.540(3) 


3.357(6) 


1.759 


0.74 


0.473 


2.2 
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0.4017(2) 


2.415(5) 


5.11(2) 


1.720 


1.17 


0.506 


1.35 
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0.4337(3) 


3.42(1) 


7.27(4) 


1.707 


1.76 


0.527 


1.00 
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0.4596(4) 


4.52(2) 


9.60(11) 


1.682 


2.90 


0.548 


0.98 
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0.5008(5) 


7.05(2) 


15.9(4) 


1.690 


4.62 


0.582 


0.62 


10 


0.5343(6) 


9.90(3) 


23.2(11) 


1.691 


7.0 


0.600 


0.50 


12 


0.5629(8) 


13.15(6) 


34.(4) 


1.70 


10.6 


0.610 


0.53 


14 


0.588(1) 


16.71(8) 
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14.1 


0.62 


f» 0.6 


16 


0.612(2) 


20.54(10) 
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19.0 


0.65 
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18 


0.632(2) 


24.73(14) 
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20 


0.652(2) 


29.3(2) 
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26. 


0.64 


^ 0.5 


24 


0.689(3) 
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39. 
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fa 0.7 
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54. 
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^ 0.7 


35 
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76. 
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0.790(4) 


94.6(20) 
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0.846(5) 
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0.870(7) 















distances. 

From general scaling arguments we expect U(r) and 
V{r) to depend on the arm length N, for ^ 1, only 
via the scaling variable x = ^/Rg, where Rg is the 
gyration radius of the star (for large N, Rg scales as 
Rg w yZ^iV, with values of A/ given in Table 1). We 
checked this by making plots similar to Fig. 1 also for 
other values of N (not shown here) and by estimating 
U{0) for different N. As argued in ja, HI for linear poly- 
mers (/ = 2), the convergence for iV — > oo is from above, 
Unj{0) « C/oo,/(0) + for small / (/ < 6); for 

larger / the data were ambiguous. For / = 2 and 4, 
which are the only cases where precise comparisons to 
previous work are possible^jhe data shown in Fig. 1 are 
in perfect agreement with 8,'9]. Values of C/(0), extrap- 
olated to iV ^ oo, are also given in Table 1. They seem 
to scale as U{0) w 0.6f^-^^. 

Much more attention had been given in the literature 
previously to the potential V{r) with r being the cen- 
tral monomer distance, and we shall also concentrate on 
V{r) in the following. The first result on it was obtained 
by Witten and Pincus 0|. They pointed out that the 
scaling 

4j - ^^'^''N'^'-' , (2) 

of the partition sum of a star with / arms and arm length 
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FIG. 2: Logarithms of partition functions Z)^ j:{r) against A'^ 
for / — 12. The data collapse expected from the cross-over 
ansatz Eq. (|3J is shown in the insert. 

N, together with the assumption that is 
for any fixed / a function of x = r / Rg only, 

4jM/[4j]'='A/(r/i?,), (3) 

implies that 

Vir)^Vwpir)^bfHafRg/r) (4) 
for 1 < r < i?g with 

^/ = (27/ -72/ - l)/t^ ■ (5) 

Precise estimates of 7/ can be found in They show 
that the scaling bf ^ f^^^ obtained in by assuming 
the phenomenological Daoud-Cotton model is not 
exact, a power law fit gives instead bf fa 0.27/^-^*. Both 
a/ and bf should be universal and should not depend on 
the specific microscopic realization. 

(2) 

This is illustrated in Fig. 2 where we show In Z)^'j{r) as 
a function of iV, for / = 12 and for three different values 
of r. In contrast to the data shown in Fig. 1, these data 
were obtained by growing the two stars at distance r and 
with the mutual interactions taken into account during 
the growth 0. This allows to measure Z^^j:{r) down to 
very small distances and large N, where it is so small that 
the ratio Z^y(r)/[Z^y-]^ measured from independently 
grown stars would be indistinguishable from zero. On the 
other hand, at large distances this second method would 
give very bad estimates of V{r), since it is obtained by 
subtracting the (nearly equal) free energies obtained in 
two independent runs. Therefore, in the following, all 
plots will show data obtained either by the first or by 
the second method, or will contain combinations of both 
types of data. 

Eq. 10} cannot hold for large distances, and it is there 
where previous results were most uncertain. An analytic 
ansatz which is supposed to cover all values of r was made 
by Likos et al. ^J. Using a 'corona' radius a which is 
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FIG. 3: V{r) for / = 18, plotted on a logarithmic scale against 
r/Rg. The short continuous curve on the left corresponds to 
Eq. Q, the dashed curve is a Gaussian. 



3.5 
3 
2.5 

N 

^ 2 



N 



1.5 

1 

0.5 




r/R„ 





f% '■=10 




-% r = 20 X 




1.+ r = 40 □ 


/ 


1 Eq.{6) 







FIG. 4: Rescaled radial Mayer functions against r/Rg for 
/ = 18. The dashed curve is the prediction from Eq. © with 

a = l.SRg 



roughly comparable in size to the gyration radius, they 
assumed that 



V{r) = 
18 



— ln(r/a) 

a j T 



-\/7/5 



1+V7/2 

exp (-Vf{r 



for r < a 



r) /2(t) for r > ct 



(6) 



This was supported by molecular dynamics simulations 
and was also shown to be compatible with experimental 
results. It was used in extensive simulations of semi- 
dilute and concentrated solutions, and gave rise to a 
number of very interesting predictions 0, 0. But for 
linear polymers it disagrees with the analytic results of 
and seems hard to be reconciled with the simula- 
tions of H, |3- In particular, it was shown in 
that Eq. © is in gross violation with simulations of off- 
lattice stars with up to 18 arms. But these arms were 
very short^whence one might doubt the relevance of the 
results of 1S| . 

Anyhow, in a later paper Jusufi et al. 19] proposed to 
use Eq. © only for / > 10, and to replace it for / < 10 
by an ansatz with Gaussian decay for r > a, 



V{r) 



A f3/2 
18 



-ln(r/cr) + 375^ 
2^exp(-T2(r2 - 



for r < a 
a'^)) for r > a 



(7) 



Notice that this does not alleviate the serious conflict 
with jl^]. Also, we would expect that the center of mass 
gets closer to the central monomer as / increases. Thus, 
if U (r) is roughly Gaussian for large r, we should expect 
that also V{r) is Gaussian there for / ^ 1. 

Let us for the moment concentrate on / = 18 arms, the 
case studied in . In order to get a first overall impres- 
sion of V{r), we show in Fig. |31its logarithm, obtained 
for fixed r = 20 and for all N < 400, against r/Rg. The 
short continuous curve at small r is the Witten-Pincus 



prediction, modified by taking the measured values of 7/ 
and 72/. It is relevant only for r <C For r ^ Rg the 
MC data can be approximated by a parabola, i.e. V{r) 
is roughly Gaussian, 



V{r) W VGa,nss{r) = Cfe 



(8) 



We conjecture that c/ and df are universal. A Yukawa 
tail as in Eq. © would essentially correspond to a 
straight fine in Fig. 13 and is definitely ruled out p^. 

Since r = 20 is not very large, one might be worried 
about finite size corrections. When plotted as in Fig. O 
finite size corrections would be visible only in the r.h.s. 
tail where V is so small that they are irrelevant. Thus 
we plotted in Fig. ^ the rescaled radial Mayer function, 

(r/RgffMir) - (r/Rgfil - exp(-y(r)) (9) 

which is the most interesting quantity, for three values 
of r. This plot agrees very well with the simulations 
of fl^, although those authors used a continuum model 
with much shorter arms. On the other hand, our data 
disagree strongly with Eq. which is indicated by the 
dashed curve. 

Linking the small- and large-r behaviours seen in Fig.|21 
into a piecewise analytic form as in Eqs. (|SJ) or 10 would 
obviously give a discontinuous slope and a bad fit. Rather 
we found that the following ansatz describes all our data 
quantitatively, for all 2 < / < 35 and for all values of r: 



V{r) = — In 



(10) 



with Vwp{r) and VGauss(?') defined in Eqs. (@J and ((HJ, 
and with r/ being an additional parameter for every 
/. It is easy to see that V{r) > for all r and that 
Vir) = 1/Gauss(r) [1 + Oir-'^f)] for r ^ 00, while 
V(r) ^ Vwp{r) [1 + 0(r2)] for r ^ 0. Like the pre- 
vious parameters, also Tf should be universal. Values for 
af,Cf,df, and r/ obtained by fitting our MC simulations 
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FIG. 5: Rescaled radial Mayer functions against r/Rg for 
several values of /. Curves are obtained from Eq. IIUII . with 
fitted parameters af,Cf,df and r/ given in Table 1. 

are given in Table 1. One sees that r/ is between 1/2 and 
1, except for the smallest values of /. The strength of 
VGauss('') increases roughly as c/ « 0.1/^-^*. Its range in- 
creases faster than Rg and the peak of the radial Mayer 
function increases even faster, roughly as Rglnf. For 
several values of /, radial Mayer functions arc shown in 
Fig. El together with the fits obtained with Eq. (|10|l . For 
/ 3> 1, their peaks are at r/Rg > af, i.e. at distances 
where Vwp(r) would be negative. 

For / > 20 our ansatz for V{r) can be simplified. 
For such stars the potential is so big for small r that 
the Witten-Pincus term can be neglected for dilute so- 



lutions: Whenever it would be relevant in comparison 
to the Gaussian term, the pair distribution function 
exp(— y(r)) is already zero for all practical purposes [2^ . 
The Witten-Pincus part becomes important only for very 
dense systems. But there the description in terms of ef- 
fective two-body forces is questionable. For the same 
reason, also the parameter rj is less precisely determined 
than af,Cf and df. 

In summary, we have obtained very precise Monte 
Carlo estimates of the effective potentials between two 
star polymers with equal number of arms and equal arm 
lengths. Using a soft core polymer model (the Domb- 
Joyce model) we have reduced corrections to scaling to 
a minimum, and we have been able to simulate many 
arms without having to use a large central particle. We 
thus believe that our results present essentially the scal- 
ing limit of long arms. Our most important finding 
is that effective potentials are much harder than previ- 
ously believed. This refers to the case where the central 
monomers are used to define the distance. For the alter- 
native case of the center-of-mass distance, it had already 
been assumed by previous authors that the potential is 
relatively hard at large r and approximately Gaussian. 
We found that basically the same is true also for the cen- 
tral mass definition. Which of these two alternatives is 
a better starting point for effective potentials in systems 
with finite concentration is another question, but our re- 
sults suggest that it does not make much difference. 

We thank Walter Nadler for discussions and for criti- 
cally reading the manuscript, and Christian von Ferber 
for interesting correspondence. 
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